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Chapter 1 



Introducing Two-Dimensional 
Bose Gas 

1.1 Low-dimensional Bose gases 
1.1.1 General overview 

The exotic phenomena of Bose-Einstein condensation was first predicted by 
Einstein , generalizing the concept of earlier work of Bose j2] on the quan- 
tum statistics of photons to indistinguishable particles with integral spin. 
These particles, later known to be bosons obey the symmetrical property of 
wave function under the exchange of particles. In nature atoms are classified 
into bosons or fermions according to their spin properties or intrinsic angular 
momentum. Bosons have integral spin while fermions are half-integral spin 
particles and are described by anti-symmetric wave function. Pauli's exclu- 
sion principles prohibiting two fermions to occupy a same quantum state is 
also due to this anti-symmetric property. On the other hand the symmetric 
nature of bosons allow them to macroscopically occupy a given state. 

The most illustrative example of bosonic behaviour of atoms is the process 
of Bose-Einstein Condensation (BEC) where at sufficiently low temperatures 
the bosons would become locked together in the lowest single-particle quan- 
tum state. Bosons that were distributed among many eigenstates plunge 
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macroscopically into the equilibrium ground state (lowest energy state) via a 
phase transtition. At this stage the de Broglie waves of neighbouring atoms 
coalesce to form a giant matter-wave. For the three dimension (3D) geome- 
try the condition for attaining BEC is to achieve a dimensionless phase space 
density of 

p\ 3 DB > 2.612 (1.1) 

where p and Xdb are the density of the system and de Broglie wave length 
respectively. One can also say that BEC occurs when the de Broglie wave- 
length becomes comparable with the interparticle distance d = p" 1 / 3 . 

BEC has been actively pursued over the last few decades in a variety of 
different systems. The initial drive in studying BEC focused on the intepre- 
tation of the collective (superfluid) behaviour of liquid 4 He. In this case 
however, the interactions between atoms are so strong that the depletion 
of the condensate is very large. To investigate the physics of BEC phase 
transition in detail experimentalists have sought to obtain the phenomena 
of BEC in a weakly-interacting systems, where the interactions do not mask 
the quantum statistical effects. 

The first "true" Bose- Einstein condensate was created by Cornell, Wieman 
and co-workers at JILA on June 5, 1995 [3J. They did this by cooling a 
dilute vapor consisting of approximately 2000 rubidium-87 atoms to 170 nK, 
the lowest temperature ever achieved at that time, using a combination of 
laser cooling (a technique that won its inventors Steven Chu, Claude Cohen- 
Tannoudji, and William D. Phillips the 1997 Nobel Prize in Physics) and 
magnetic evaporative cooling. The Laser cooling was followed by cycles of 
evaporative cooling, a principle similar to the thermalization of a cup of 
coffee. Henceafter a peak appears at the centre of the density distribution of 
the atomic cloud, which can be imaged by a resonant absorption technique 
after a ballistic expansion of the cloud as it is realesed from the trap (See 

Fig. (HID). 

About four months later, an independent effort led by Wolfgang Ketterle at 
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MIT created a condensate made of sodium-23 jl] . Ketterle's condensate had 
about a hundred times more atoms, allowing him to obtain several important 
results such as the observation of quantum mechanical interference between 
two different condensates. Cornell, Wieman and Ketterle won the 2001 Nobel 
Prize for their achievement. In that same year Bose condensate on lithium 
was also observed by the group at Rice led by Randy Hulet [H] 

The initial results by the JILA, MIT and Rice groups have led to an explo- 
sion of experimental activity. For instance, the first molecular Bose-Einstein 
condensates were created in 2003 by teams surrounding Rudolf Grimm jB] at 
the University of Innsbruck and the first fermionic condensate by Deborah 
Jin and co-workers [Zj at the University of Colorado at Boulder. 

These pioneering realisation of BEC in a weakly-interacting system of 
trapped dilute alkali atoms in the mid 90 's has revived the interest in the 
theoretical studies of Bose gases. A rather massive amount of work has been 
done in the last couple of years, both to interpret the initial observations 
and to predict new phenomena. In the presence of harmonic confinement, 
the many-body theory of interacting Bose gases gives rise to several unex- 
pected features. This opens new theoretical perspectives in this interdis- 
ciplinary field, where useful concepts coming from different areas of physics 
(atomic physics, quantum optics, statistical mechanics and condensed-matter 
physics) are merging together. 

Due to the enormous breadth of the subject, it is quite impossible to give 
a full coverage of the basic theory of Bose-Einstein condensation and its 
experimental realizations together with its many novel applications. A large 
series of reviews and books has already emerged elaborating in clear details 
this novel state of matter. Readers are recommended to retrieve informations 
from Pethick and Smith jHj, Legget [Oj, Parkin and Walls Griffin, Snoke 
and Stringari and so on. This thesis is written for a specialist audience 
for those whom already acquinted with the basic theory of Bose-Einstein 
condensation in a 3D geometry. In spite of this, or rather because of this I 
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Figure 1.1: The first BEC in a gas of 87 Rb atoms, as it appears in two- 
dimensional time-of-flight images as temparature is decreased below T c . 
Adapted from Anderson et. al. [3] 

move ahead bypassing into the two-dimensional (2D) territories introducing 
the peculiarity phenomena arising in the squeezzed flat geometry. 

1.1.2 Bose-Einstein condensate in the flat geometry 

The possibility of Bose-Einstein Condensation (BEC) phase transition occur- 
ing in a 2D geometry has a remarkable history and breakthroughs. The effect 
of dimensionality on the existence and character of BEC or superfluid phase 
transition has seen a spurt of interest in recent years. A peculiar feature 
of the low-dimensionality (ID or 2D geometry) at finite temperature is the 
absence of a 'true condensate' following the Bogoliubov k~ 2 theory, believed 
to be originated from the long wave-length phase fluctuations ^2] • This has 
been proved in the seminal work of Mermin and Wagner ^3] and Hohenberg 

mi 
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Generally for a weakly interacting Bose gas at rest, the phenomenon of 
BEC is related to the macroscopic population of the single-particle state at 
zero momentum. One should expect a peak in the momentum distribution 
if condensation occur. One can also rely on the criterion first suggested 
by Onsager and Penrose ^Hj in the context of quantum field theory. They 
associated BEC with the emergence of long-range "off-diagonal" order in the 
one-body density matrix p(r — r') =< i/)<(r)ip(r') >. Here ijj(r) and ifj '(r) 
are the annihilation and creation field operators for a particle at position r 
and < ... > is the mean value. It was shown by Hohenberg [Hj that for two- 
dimensional Bose gas the correlation function depends algebraically on |r— r'| 
with an exponent 7(T) that varies proportionally with the temperature T as 



Algebraic decay of the correlation function is what one expect when the 
temperature is tuned to the critical temperature of a continuous phase tran- 
sition. In fact, this is known as the Kosterlitz-Thouless transition where the 
system cools toward a superfluid state through the binding of vortices with 
opposite vorticity ^H]. Nevertheless in the limit T — > global coherence is 
achieved in a homogeneous 2D system and a true condensate then exists. 
In contrast, for a trapped 2D fluid the modification of the density of states 
caused by the confining potential allows a true condensate to exists even at 
finite temperature. 

Recent experiments have realized a quasi-two-dimensional (Q2D) cold Bose 
gas by tuning the anisotropy of the trapping potential [T71 UHl EH 120] • This 
achievement has attracted vast attention on the new physics appearing in the 
currently studied low- dimensional trapped gases and the nature of finite-size 
and finite-temperature effects. 





5 



Low- dimensional Bose gases 



1.1.3 Homogeneous non- interacting Bose gas 

For a homogeneous gas, the individual particle eigenstates are characterized 
by momentum k and corresponding eigenenergy ek = h 2 k. 2 /2m, where m 
is the mass of the particle. At thermal equilibrium with temperature T, a 
uniform gas of N bosons obeys the well-known Bose-Einstein distribution 

/(ek, T) = 1 , (1.3) 

exp((e k - n)/k B T) - 1 

where k B is the Boltzmann constant and \x corresponds to the chemical po- 
tential of the system. The chemical potential can be interpreted as the energy 
required by a N-particle system to exchange one particle with its environ- 
ment. The total number of particles in the excited states is given by the sum 

iVex = ^/(e k ,T), (1.4) 

k 

which specifies the value of the chemical potential \i. However, in a macro- 
scopic system the volume is so large that the energy levels can be treated as 
continuous. Thus, Eq. (|1.4jl can be written in the integral form 

N ex = r j de k D(e k )/(e k ,T) (1.5) 

with r being the two-dimensional area and -D(e k ) is the density of state. In 
two-dimension, the smoothed density of state is D(e k ) = L 2 m/2Tch 2 , where 
L 2 is the surface area. If we take the lower limit to be zero, then the integral 
does not exist because of the divergence at low energies. In a two-dimensional 
box, the ground state has energy eo ~ h 2 /mL 2 , since its wavelength is com- 
parable to the linear extent of the area. To estimate the number of excited 
particles, we implement a cut-off (eo) in the integral as a lower bound. This 
gives us 

k B T k B T 

N ex ln( ). (1.6) 

The number of excited particles becomes equal to the total number of parti- 
cles at a temperature (critical temperature) given by 

kBTc ~^\^N =h ^N> (L7) 
6 
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where a = N/L 2 is the number of particles per unit area. For a large sys- 
tem and in the thermodynamic limit (N —* oo), the transition temperature 
tends to zero. This indicates that the macroscopic occupation of particles 
can only occur at zero temperature for a uniformly large system in the ther- 
modynamic limit. Hohenberg pi] has proposed an ad absurdum reasoning 
for the absence of off-diagonal long range order for in a 2D uniform Bose gas 
at finite temperature. In the presence of condensate the following inequality 

n w 1 mk B Tn 

n(k) ^-2 + W7T (L8) 

should hold, where n(k) is the momentum distribution at low momenta, Uq 
is the condensate density while n is the total density. However for 2D, the 
density of excited particles evaluated using the integral n exc = f d 2 kn(k) 
diverges. So, one is led to the conclusion that the condensate fraction must 
be zero at any finite temperature. 



1.1.4 Kosterlitz-Thouless transition 

The theorem of Mermin and Wagner ^3], which forbids the spontaneous 
breaking of a continuous symmetry in two-dimension due to the enhanced 
importance of long-wavelength thermal fluctuations, also reinforces the argu- 
ment put forward by Hohenberg [T^j. However, that a different kind of phase 
transition to a state with algebraic long-range order occurs in a 2D Bose 
fluid was first pointed out by Kosterlitz and Thouless jTHj. This so called 
topological phase transition is associated with dissociation of a large number 
of vortex pairs and the immediate destruction of superfluidity caused by the 
phase slip process that occur once unbound vortices are present [23 • Nelson 
and Kosterlitz [22] predicted by means of a Renormalization-Group analysis 
that the discontinuous drop in the superfluid density n s at the critical tem- 
perature T c is of universal nature and such that n s A 2 = 4, where A c is the 
thermal de Broglie wavelength. 

The Kosterlitz-Thouless transition can easily be demonstrated using the 



7 



Low- dimensional Bose gases 



two-dimensional XY-model. The model consists of planar rotors of unit 
length arranged on a two-dimensional square lattice. The Hamiltonian of 
the system is given by 

h = -j s * • s j = - J cos (^ ~ fa) ( L9 ) 

<i,j> <i,j> 

where J > and the sum < i,j > is over the lattice sites of nearest neigh- 
bours only. Here we take |Sj| = 1 and 0; is the angle that the ith spin makes 
with some arbitrary neighbours. Expanding about a local minima of H gives 

H-E ~iJ dr(V0(r)) 2 (1.10) 

where 0(r) is a phase function defined over the lattice sites and E is the 
non-vortex ground state. By spherical symmetric consideration and taking 
into account a quantized vortex configuration, we can estimate the energy of 
a single vortex E vor to be 

E vor -E ~7rJln(L/a s ) (1.11) 

where L is the finite system size and a s is the lattice spacing. However, 
the free energy of a single vortex can be calculated using the Helmholtz free 
energy relation F = E vor — TS, where S is the entropy estimated from the 
number of places where we can position the vortex center (in a classical 
viewpoint), namely on each of L 2 plaquette in a square lattice, i.e. S = 
ln(L 2 /a 2 ). Accordingly the free energy is given by 

F = E + (vrJ -2k B T)\n(L/a s ) . (1.12) 

For T < 7r J/2ks the free energy will diverge to plus infinity as L — > oo. 
At temperatures T > n J/2ks the system can lower its free energy by pro- 
ducing vortices (F — > — oo). This simple heuristic argument points to the 
fact that the logarithmic dependence on system size of the energy of the 
vortex combines with the logarithmic dependence of the entropy to produce 
the subtleties of the vortex unbinding transition. It is the logarithmic size 
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dependence of the two-dimensional vortex energy that allows the outcome of 
the competition between the entropy and the energy to change qualitatively 
at a certain finite critical temperature T^t known as the Kosterlitz-Thouless 
transition. 

In reality it is not the single vortices of the same signs that proliferate at 
a certain temperature. What happen is that the larger vortex pairs which 
are bound together for temperature below Tkt unbind at Tkt- This is a 
collective effect. 



1.1.5 Non-interacting trapped Bose gas 

Let us consider a symmetric harmonic confining potential V(r) = muj 2 (x 2 + 
y 2 )/2 to describe the BEC cross-over of the 2D Bose gas. In this case the 
particle energy is E v = h\oin x + n y ), with quantum numbers n x and n y 
which are non-negative integers. The density of states for this system is then 
D(E) = E/(ftw) 2 , and the contribution of low-energy excited states to the 
total number of particles is negligible. The ground state (E = 0) population 
is given by 

exp(-n/k B T) - 1 

which becomes macroscopic even for a small negative \i. A cross-over to BEC 
regime exists in this case. Therefore, separating out the population of the 
ground state, one can replace the summation in Eq. (jl.4j) by an integration, 
and one gets 

N = N + V J rfe k D(e k )/(e k , T) . (1.14) 

Using a the relation —{i/hsT w 1/iVo, which is obtained using Eq. ()1.13j) 
assuming that ground state is macroscopic populated one can easily show 
that Eq ()1.14j) gives the following result 



N 



where T c = ^J^Huj is the 2D BEC transition temperature for the non- 
interacting case. The effects of the interactions on this critical temperature 
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will be discussed in chapter 2. 

1.2 Weakly interacting Bose gas 

Preliminary works of Bogoliubov [21], Huang et al. [213 121] and Lee et al. 
|2T[ |2H] were mainly focused on the theory of interacting homogeneous Bose 
gases. The extension to inhomogeneous gases was first explored by Gross j2H] , 
Pitaevskii [SU] an d Fetter p?T) l32*]. The thermodynamic behaviour of trapped 
Bose gases has been already the object of several theoretical works, starting 
from the development of the Hartree-Fock formalism [23], the inclusion of 
interaction effects using the local density approximation p^H IHH] up to the 
most recent approaches based on self-consistent mean-field theory [TTJ EE] 
and numerical simulation. 

1.2.1 Weakly interacting regime 

I will consider weakly interacting gases with a short-range potential of in- 
teraction between particles. The total interaction energy is equal to the 
sum of pair interactions in the c-number formalism and can be written as 
Em = N 2 g/2A, where g is the coupling constant for the pair interparticle 
interaction, N is the number of particles, and A is the area. Consequently, 
the interaction energy per particle is equal to E int = ng, with n being the 2D 
density. The general criteria for a weakly interacting regime assumes that 
the mean interparticle separation f greatly exceeds the characteristic radius 
of interaction between particles, l e . In the 2D case we have f ~ (27m)" 1 / 2 
and thus one obtain the inequality nil ^ 1- 

Another important condition for weakly interacting regimes is that the 
wavefunction is not influenced by the interaction between the particles if 
the interparticle distance is of order f. Based on this condition, we can 
develop a physical picture which will be used for finding how the criterion of 
the weakly interacting regime depends on the dimensionality of the system. 
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For example let consider a box of size f, which accomodates one particle in 
average. In the limit T — > 0, the particle kinetic energy is Ekin ~ h 2 /mf 2 . 
The wavefunction is not influenced by the interparticle interaction if the 
kinetic energy per particle is much larger than the interaction energy per 
particle. The inequality E^ n ^> Ei nt /N immeadiately gives the criterion of 
weakly interacting regime in term of of the density and coupling constant. 
In the 2D case this criterion takes the form 



where g is the 2D coupling constant which will be explained in detail in 
the next chapter. Most of the work in this thesis deals within the small 
repulsive/attractive limit (g > or g < 0), except for the case where a 
weakly bound state of colliding atoms is present near resonance. Then the 
strength of g is extremely large and even at very low densities the criterion 
of Eq. (jl.l6|) does not hold. Thus, one has to look for solutions beyond 
mean-field. 

1.2.2 Elementary excitations and quasiparticles 

The presence of a condensate mean field is naturally expected to affect the 
spectrum of the elementary excitations of a bosonic gas. This is a direct con- 
sequence of the fact that the atoms are not interacting with other individual 
atoms, but they are moving dominantly in the presence of a coherent con- 
densate field. As expected, the deviation from the non-interacting spectrum 
become more pronounced as more particles are added to the condensate, 
since this lead to an increase in density. The notion of an excited (bare) 
atom interacting with condensed atoms leads quite generally to the idea of a 
dressed atom, or quasi-particle. One can see the significance of interactions 
to the change in the spectrum of elementary excitations by looking at the 
quasi-particle dispersion relation given by [21] 



m\g 



< 1. 
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where = h 2 k 2 /2m corresponds to the non-interacting energy spectrum, n c 
is the condensate density and Vq is the Fourier transform of the interaction 
potential in the limit k — > 0. 

Recently, it was demonstrated by Al Khawaja et. al [HZ| that by taking 
the following modifed type excitation spectrum: 

fuj k = yje 2 k + 2fie k (1.18) 

where \x = n c T 2B (—2[i), one is able to avoid the dangerous ultraviolet and 
infrared divergences in the 2D Bose gas. T 2B (— 2/i) is the multi-loop two- 
body T-matrix taken at energy argument —2/i which is exactly what it costs 
to excite two atoms out of the condensate (See Sec. 1.3 for more details). 



1.2.3 Macroscopic framework of weakly interacting trapped 
Bose gas 

Let us consider the Bose gas trapped in a 2D isotropic planar trap V(r) = 
muj\{x 2 + y 2 )/2. The starting point is the grand-canonical Hamiltonian in 
second quantization language: 

H = H + ^JdrJ dr'^(r)^(r')V(r - r')^( r ')^( r ) ( L19 ) 

where H = J dr^(r) — |^V 2 + V(r) — /x ip(r) , /x is the chemical poten- 
tial and V(r — r') is the interatomic interaction potential. The annihilation 
and creation operators are denoted by ip^(r) and ip(r') respectively and obey 
the Bose- Einstein commutation relation 

[^(r),^ f (r')] = 5(r-r') [^(r), ^(r')] = [^(r), ^(r)] = . (1.20) 

Since we are looking at system at very low temperature, only the contribution 
of s-wave scattering is important. Consequently we neglect the momentum 
dependence of the interatomic interaction and use V(r — r') = Vq5(t — r'). 
This leads to ultraviolet divergences but can be easily remedied. I will illus- 
trate the technique in the next chapter. However with this contact potential 
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the above Hamiltonian will reduce to the following one, 

H = H + ^J drV^M^M^r^r). (1.21) 

The time- dependent Heisenberg operator ip(r, t) = exp(iHt / h)ijj(r) exp(—iHt/h) 
obeys the equation of motion 

ih^^-^mr,t),H] (1.22) 

which yields the non-linear operator equation 

= + V(r)j 0(r, t) + g^(r, t)^(v, t)^(v, t) . (1.23) 

By separating out the condensate part of the particle field operator one can 
write ip(r) as a sum of spatially varying condensate wavefunction \l/(r) and 
a fluctuation field operator z/>, 

= $(r)+^. (1.24) 

The wave function \l/(r) is also known as the order parameter and is defined 
as the statistical average of particle field operator \l/(r) =< ip >. 

Using the decomposition in Eq. (jl.24)l above, one can write the interaction 
term ip\r, t)ip(r, t)ip(r, t) of Eq. (jl.23j) in the form 

f (r, t)ip(r, t)tp(r', t) = |^| 2 ^ + 2|^| 2 ^ + ty 2 ft + + 2^ip + . 

(1.25) 

Taking the average of the above relation one obtains 

< ^(r, t)ip(r, t)^(r, t) >= + m^* + 2n nc + < ip ] ipip > , (1.26) 

where n c = |\1/| 2 is the condensate density, n nc =< ip^ip > is the non- 
condensate density while rh =< ipip > is the off-diagonal (anomalous) den- 
sity. So taking the average for Eq. (jl.23j) and using the relation in Eq. (jl.26|) 
we arrive at the equation of motion for \l/(r,t): 

d^(r t) ( h 2 \ 
ih m = l-—V 2 + V{r)+gn c {r,t) + 2gn nc {r,t))y{r,t) 

+ grh^* + g < ^ ] 4>4> > (1.27) 
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Neglecting the three-field correlation function < ip^vpip > in the equation 
above is called the Hartree-Fock-Bogoliubov (HFB) approximation. Further 
simplification to that by ignoring the term m =< ipip > is known as Hartree- 
Bogoliubov-Popov (HBF-P) approximation. In this thesis I will refer to the 
static Popov limit in which the fluctuation of the thermal cloud is ignored 
n nc (r,t) ~ n nc (r) and will only consider a stationary solution for \I/(r,i) = 
^(r) exp(— ifit/h). Thus the Non-linear Schrodinger equation Eq. (jl.27j) 
leads to the Gross-Pitaevskii (GP) equation 

h 2 \ 
-—V 2 + V(r) + g 2 n c (r) + 2 9l n nc (r) J tt(r) = /^(r) . (1.28) 

As anticipated, the general finite temperature GP equation above is not 
closed. Its solution requires one to know the information about the distribu- 
tion of the non-condensate particles. Readers should also be aware that in 
the above equation I have introduced g\ and g<i which represent the coupling 
parameters due to condensate-noncondensate and condensate-condensate re- 
pulsion respectively. This choice of coupling parameters in studying the 
equilibrium property of condensate-thermal particles is the core issue of this 
thesis. Coming back to the problem of thermal density distribution, I will 
take a semi-classical approximation [2H] supplemented with a Hartree-Fock 
type energy spectrum 

v 2 

e(p, r) = f - + V ext {r) + 2 9l n c {r) + 2 9l n nc r) . (1.29) 
Hence, the thermal density is obtained by integrating out the momentum 

nUr) = J (2^ /(e(p ' r) ' T) 

' xln[l-exp03Gu-K//(r))] . (1.30) 



27Th 2 

where we consider the effective potential as V e //(r) = V(r) + 2gin c (r) + 
2gi n nc (r). This type of method that couples solution of condensate and 
non-condensate density is also known as the Two-Fluid model. 
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1.3 Brief introduction of the T-Matrix 



To describe the microscopic treatment of binary collisions, one generally 
considers the interaction between two condensed atom in a dilute weakly- 
interacting BEC. The complete set of two-particle interactions can be summed 
up into a single operator called the Transition matrix or simply T-matrix. 
This operator can be determined via a perturbation theory expansion in term 
of the actual interatomic potential. The possible scattering events can be de- 
scribed by a single interaction (k'| V |k), or alternatively the particles may 
first make a transition to an intermediate state |q) before interacting again 
to emerge in the state |k'). The repeated form of the scattering from initial 
to a final state via the intermediate state |q) can be defined as an expansion 
of the T operators in term of the actual interatomic potential V via 



where G° = (z — H°)~ l corresponds to the unperturbed Green's function and 
H° is the free Hamiltonian in the absence of interactions while z corresponds 
to the centre-of-mass-energy of for the pair of atoms. In practice, the T 
operator can be assembled more explicitly in term of the Lippman Schwinger 
relation (see also in Sec. 2.1.1). 



The two-body interaction above describe collision in vacuo in which the 
intermediate states are single-particle in nature. In order to illustrate the 
interaction within a two-dimensional BEC we must consider the scattering 
in the presence of condensate particles. Interactions of this nature can be 
described by the many-body T-matrix T MB (k, k', K; z) which due to the 
surrounding particles depends also on the center-of-mass momentum K of 
the colliding particles. In the zero momentum limit we set k', k, K = 0. 
Since we are interested in the matrix element with all momenta equal to zero 
we find to a good approximation the following Bethe-Salpeter equation for 



T = V + VG°V + VG°VG°V + .... 




T = V + VG"T . 
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the scattering of a quasi-particle, 



T MB (k,0,0;z) = V(k) + -XV(k') 



k' 



1 + 2N(nw kf ) 

Z — IfrhJ^j 



x T MB (k',0,0;^) 



(1.33) 



where hv-^ is the quasi-particle energy defined in Eq. (jl. 18)1 and N(x) = 
l/(exp(/?x) — 1) is the Bose-Einstein distribution function. However in the 
zero momentum limit, the above Bethe-Salpeter equation is tainted by an 
ultraviolet divergence due to the neglect of the momentum dependence of the 
potential V(k). This divergence is remedied by using the following relation, 



where z = — j/i with j = 1,2 depending whether the collision is between two 
condensate atoms (j = 2) or between one condensate and one non-condensate 
atom (j = 1). The details are given in Sec. 2.1. The collision between 
condensate and the thermal atoms become more probable as temperature 
increases. It is a crucial point and covers our investigation on the effect 
of interaction in calculating the density profiles at finite temperature (See 
Chapter 2). 

1.3.1 Quantum control of the interacting Bose gas 

A Feshbach resonance occurs when the energy of a metastable molecular state 
E res matches the energy of a pair of incoming condensate atoms E t h and a 
coupling exists between them during a collisional process [HE]- As described 
in the previous section, the scattering of particles with zero relative energy or 
momentum is specified by the scattering amplitude related to the T matrix. If 
the energy of the incoming particles Eth is close to the resonance energy E res 
of one particular bound state (resonance state \ip res )), the contributions from 
all other states will vary slowly with energy and they may be represented by 
non-resonant scattering length a^ g whose energy dependence is much weaker 




(1.34) 
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compared to the resonance term. By writing the T matrix elements as a 
function of scattering length and using the argument discussed above we 
obtain the following relation, 

T MB (0, 0, 0) = T b A / B (0, 0, 0) + C{B) . (1.35) 

&th ~ J^res 

where T^ B is the off-resonant scattering amplitude and C (B) is the resonant 
width, weakly dependent on B. The total scattering amplitude that includes 
resonant and non-resonant terms for a 3D Bose can be easily calculated 
as T MB (0,0,0) ~ T 2S (0,0,0) = Aiitfa/m by the two-body scattering of 
particles within the Born approximation |Hj. Eq. (J1.35)) displays the energy 
dependence of the two-body T matrix and is characteristic of a Feshbach 
resonance. Hence the atomic interactions may be tuned by exploiting the 
fact that energies of the states depend on external parameters among which 
are the strength of a magnetic field. Expanding the energy denominator 
about the value of magnetic field we find 

E th - E res « (fi res - 2/i i )(£? - B ), (1.36) 

where /ij = — |^,z = l,2 are the magnetic moments of the two atoms in the 
open channel while /i res = — dEr ^ s is the magnetic moment of the molecular 
bound state. Substituting Eq. (jl.36|) into Eq. (jl.35j) one obtains [21] 

a(B) = a bg (l - ^5L) . (1.37) 

where a\, g is the off-resonant scattering length. By sweeping across the mag- 
netic field around the resonance location B the coupling strength of the con- 
densate can be used as a controlling tool. Hence one has the possibility not 
only to control the scattering length from a strongly repulsive to a strongly 
attractive regimes but also in obtaining a novel molecular condensate [Hj (as 
depicted in Fig (jl.2|) . 

Whether we can have a similar control of the interaction parameters for 
the strictly 2D condensate when the s-wave scattering length becomes larger 
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' " + 

Magnetic field 

Figure 1.2: Energy diagram for the atomic scattering state and the molecular 
bound state. The Feshbach resonance condition occurs near 20 G, where 
the Zeeman energy of the atomic scattering state becomes equal to that 
of a molecular bound state because of the difference in magnetic moments. 
Molecules at (2) are created from the BEC at (1) by a downward sweep of the 
magnetic field across the resonance. For detection, a reversed sweep brings 
the molecules above the dissociation limit. The inset schematicallyshows the 
molecular potential that corresponds to the open channel (lower curve) and 
the molecular potential that supports the bound state (upper curve). U, 
potential energy; r, interatomic distance. From Herbig et. al. [Oj 
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than the width of the axial trapping ^)\ IT2] will be addressed in the second 
part of chapter 2. To my knowledge, I have not encountered any attempts 
to study the Feshbach resoances for a strictly 2D Bose gas. This motivated 
me to look into such phenomena in those low- dimensional regime by using 
the similar technique conventionally used in studying them in the 3D case. 

Once a tunable coupling parameter is calculated it will be incorporated 
into the 2D dynamical GP equation to study the effect of the interactions 
on the density profiles at absolute zero temperature. I will perform a similar 
analysis as done by Kagan et.al [H] by modifying the coupling strength as a 
dimensionless physical parameter rj which is defined as the ratio between the 
mean-field interaction energy and the harmonic-oscillator level spacing. The 
interaction is of repulsive nature if 77 > 0, else it is attractive for 77 < 0. The 
magnitude 77 <C 1 indicates that we are in the weakly repulsive (attractive) 
regime and 77 ^> 1 means a strongly repulsive (attractive) for 77 > (77 < 0). 

1.3.2 The attractive interaction limit 

For an attractive interaction between particles ( 77 < ), the picture changes 
drastically. A Bose condensate with \r]\ ^> 1, in which the structure of the 
trap levels is not important, cannot be formed at all, since in this case the 
accumulation of particles in one quantum state would be associated with 
increase in energy |41j . In this case the dominant elastic pair collision will 
prevent the formation of a Bose condensate with densities {n c /n )r] ^> 1), 
since the structure of trap levels is essentially smeared out by interatomic 
interaction and there is no gap for the excitation of particles from the con- 
densate. 

As a simple illustration, by taking the TF approximation to Eq. (jl.28|) 
(assuming n nc = at T = ) in this limit (\r]\ ^> 1) and substituting 
77 = — k 2 where k is any positive number, we obtain the relation 




(1.38) 
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What we have obtained is a monotonically increasing function which is not 
physical. The use of Eq. (jl.28|) is inappropriate to analyse the behaviour of a 
condensate with attractive interaction. Consequently, we need to incorporate 
the physics of atom losses due to three-body recombination and formation of 
metastability condensate in the attractive regime as investigated by several 
authors 0U 12 ESI EH for the 3D case. 

In order to model atom loss due to three-body recombination we add a phe- 
nomenological term proportional to the density squared ^(x)) 2 with coeffi- 
cient K3/2 jlHj. The GP equation which incorporates a three-body loss term 
reads 



ih 



dV(r, t) 



H V72 . T/ , ihK 3 2 

-— V +V ext + gn c -—n c 



*(r,t), (1-39) 



dt 

where the rate coefficient K s is related to the atom loss rate given by : 

dN 



-Ksjntdr. (1.40) 

A meta-stable solution of the above Eq. (jl.39|) exist only if one chooses a 
number N such that N < N c where N c is the critical number as reported for 
the 3D case [47 ^ l4T | l4*H | I49j . Beyond this number, the system would begins to 
collapse and lead to infinite density fluctuation to form hydrodynamic 'Black 
hole' jlH]- The calculation of this critical number for the 2D case is reported 
in chapter 2. 



1.4 Dynamics of vortices in a 2D Bose fluid 

Under appropriate stabilization conditions, such as steady applied rotation, 
vortices can form a regular array. In a rotating uniform superfluid, quantized 
vortex lines parallel to the axis of rotation form a Wigner-type crystal lattice 
|50| . At non-zero temperature, dissipative mutual friction from the normal 
component ensures that the array rotates with a same angular velocity as 
the container. A triangular vortex array is energetically favoured for vortices 
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near the rotation axis of rapidly rotating vessels of superfluid helium, (see 
Tkachenko |EEJE2]). 

However, the properties of vortices in a Bose-Einstein atomic condensate 
atoms have proven to be more tractable compared to the superfluid helium. 
Following the creation of single vortex lines in an atomic condensate by a 
phase imprinting technique [S3], large arrays have been created via mechan- 
ical rotation and stirring of the condensate with laser beam [S3J ESI EE]- 
The recent observation of Tkachenko modes (collective excitations) in the 
Bose-condensed 87 Rb condensate by Coddington et al. jSZj stimulated more 
work on this direction. Congruent to that, there has been an enormous 
number of theoretical studies in analysing the properties of vortex lattices 
[SB1 EE1 EH1 EI], Tkachenco lattice oscillations (waves) |B21 EH], the transi- 
tion into a novel properties of Fractional Quantum Hall Effects (FQHE) and 
strongly correlated uncondensed liquid [HH ES] • 

Atomic Bose condensates allow one to study superfluids over a range of 
rotational frequencies Q |55] 156] 166] . Baym [H2] has classified four distinct 
physical regimes as Q increases: (1) The 'stiff' Thomas- Fermi regime, where 
Q is small compared with the lowest compressional frequencies, sk , where 
s is the sound velocity, k$ ~ 1/R is the lowest wave number in the finite 
geometry, and R is the size of the system transverse to the rotation axis. 
The system is an incompressible fluid in this regime responding effectively 
to rotation, (ii) In the limit sko <C <C ms 2 , the system is said to be 
in a 'soft' Thomas-Fermi regime. Here the vortex core size is much smaller 
than inter- vortex spacing, (iii) When Q ^> ms 2 , we enter a fast rotation 
regime where one is ensured that for a large vortex density, the intervortex 
separation is greater than the coherence length. This idea follows from the 
work of Baym jHlj who showed shrinking of vortex cores when such limit is 
reached. In this rapidly rotating regime the system enters a effectively two- 
dimensional strong field. Viefers et al. jBB] have shown that as the angular 
momentum of N particle system gets larger, particles are spread out in space 
and the 'Yrast line' (the lowest possible interaction energy for a given angular 
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momentum L) decreases. As the Yrast line decreases, the density of particle 
also decreases in parallel. Ho [HE) on the other hand has predicted that in 
this limit particles should condense into the Lowest Landau Level (LLL) in 
the Coriolis force, similarly to charged particles in the quantum Hall regime. 
At relatively low vortex density, the vortices form a triangular lattice [HIl E2J 
which is energetically favourable. There is a gap in the energy levels ~ 2co± 
when LLL acheived. The single particle wave-function in the LLL has the 
form 0j(r) ~ z l exp(— r 2 /2d 2 ± ) where z = x + iy, i=0,l,2... and d± is the 
transverse oscillator length given by ^h/muj^_. The LLL condensate wave- 
function is a linear superposition of such states : 

N 

$LLL(r) ~ Ci z ' exp(-r 2 /24) ~ ]J(z - Zi ) exp(-r 2 /24) (1.41) 

i i=l 

(iv) Eventually the vortex lattice melts [001 IHHj arid the systems enters a 
strongly correlated vortex liquid [HH EE1 EOl [HI]- The melting is somehow 
controlled by a dimensionless Lindemann parameter 

Ci = -=(r) 2 = - 

up \l v J p 

where d = p -1 / 2 is the inter-particle distance, l v = ^tt/VL is the mean inter- 
vortex spacing and n v ( n v = Q/n ) is the density of vortices. 

This parameter controls the relative strength of the quantum fluctuations 
and the melting transition into quantum Hall regime at zero temperature. 
The melting filling factor calculated using various methods [001 103 EH] is 
found to be cl = (cx,) m ~ 0.1 . The vortices are in the lattice phase if 
cl < (ci)mj otherwise in molten phase if cl > (c^) m . One could also say 
that the vortex crystal starts melting as the vortices begin to depart from 
thier lattice positions. 

However in this work, I am taking a completely diverse approach from the 
one discussed above. The difference is two-fold, firstly I am already in a 
2D geometry so that an ultra-fast rotation to achieve this geometry is not 
required. Secondly I assume my vortices are in the liquid phase. Freezing 
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and strongly correlations need to be analyzed hence after. The following 
sections will shed more light on my approach to the N vortex problem in the 
confined 2D geometry. 

1.4.1 Vortices as charged bosons 

Let us consider a system of N bosons in a condensate confined within a 
cylindrical potential (with radial u± and axial trap frequencies Q z ) and put 
into rotation with angular velocity Q. The Hamiltonian of the system reads 

HUES! 

N p 2 1 N 
Hn = Yl 2^ + 2 m ^ 2r * ~ n & ' Yi X Pi + S V ^ ~ r i) 

i=l i<j=l 

= E (Pt " m 2 Q J x ri? + \m[{^l - ny\ + ^ 
i=i 

i<j=i 

This system is identical to a system of charge-Q particles interacting under 
the influence of a magnetic field 

B* = Vx (mVLz x r/Q) = {2mVt/Q)z. 

At sufficiently low temperatures and in the fast rotation limit Q ^> ms 2 , this 
system enters an effectively two-dimensional strong field and the dissipation 
due to scattering of thermal excitations on the vortices is negligible, thus 
the vortex dynamics is conservative. The situation considered then corre- 
sponds to the Magnus-force dominated limit of vortex motion in a very clean 
superconductor |7T| 172^. 

An indirect connection to the above argument can be related to the work 
of Nelson et al. [73] who have shown that the statistical mechanics of the 
flux-line lattice (FLL) of high-T c superconductors can be studied through an 
appropriate mapping onto the 2D Yukawa boson system. Following this mag- 
nificent idea, Magro and Ceperley [TjJ [75] have performed Diffusion Monte 
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Figure 1.3: Vortices when rotation Vt — > u±. Centrifugal anf trapping energy 
compensate and spatial extent becomes large. Adapted from Bretin et. al. 

m 

Carlo (DMC) and Variational Monte Carlo (VMC) numerical calculations to 
calibrate the vortex properties for a homogeneous 2D Yukawa boson system. 
Taking a similar step, I will treat the system of N v vortices as N bosons 
confined in a 2D harmonically planar trap ^ma^r 2 performing a Density 
functional theory (DFT) calculation (See Chaprter 3.) 



1.4.2 Quantum Monte Carlo Techniques (VMC and 
DMC) 

My DFT calculation heavily relies on the accurate data of VMC calculation 
of Magro and Ceperley [23] . So I feel it is neccessary to display here some 
idea of what the method is all about. VMC method is a technique developed 
to minimize the ground-state energy of a many-body trial function. The 
energy functional E = H \ 4>t), where <pT is a family of trial functions is 
minimized by one of its particular choice of wave-function 0o known as the 
ground state. In other word we introduce a family of trial function <j)(x,p), 
depending on the parameter p and compute the corresponding energy, 

= / 0*(p,r)#0(p,r) 

and look for the choice of p = p* such that E* is the minimum energy within 
the parameter space. Obviously the success of VMC depends largely on a 
good choice of trial function. The most commonly used trial function when 
one is interested in identifying the solid or liquid phases of a system is the 



24 



Dynamics of vortices in a 2D Bose fluid 



Jastrow pair product which reads [Hj 



<Pt = JJexp[-u(r y )] JJexp[-c(ri - Z;) 2 ] 



(1.44) 



i<j i 



where rj is the position of ith of N particles, Zj are the lattice sites, 
is an interparticle distance, and tt(ry) is the Jastrow pair function chosen 
to approximate the solution of two-body Schrodinger equation. The lattice 
is assumed to be triangular since this minimizes the potential energy [77] . 
The Jastrow function contains parameters which along with c are varied to 
minimize E(p). For the liquid phases c = but for a solid c takes a finite 
value treated as an adjustable parameter. Once a right choice of variational 
function is obtained, the ground state energy is further calculated using the 
DMC method. Starting with the imaginary time Schodinger equation, 



and by introducing f(R,T) = 4>t{R)'iP(R, t) as the mixed distribution func- 
tion to the equation above we can obtain the diffusion-reaction equation , 



= A V 2 / - 2A 2 V * • U*i<fr) - & ~ E T )f (1.46) 



where A is the dimensionless DeBoer parameter defined as A = h 2 /2ma 2 e in 
which a and e are the reduced unit length and energy. El = Hc^t/^t found 
in the last term represents the upper bound on the ground state energy. The 
value of Et that gives a steady poulation is the exact ground state energy. 
Other properties such as the pair distribution functions are also calculated 
from the mixed distribution, ^t(R)^o(R) where ^o(-R) being the ground 
state in a large time evolution. The pair-correlation function g(r) and the 
structure function S(k) encode the information on the degree of order of the 
system. Oscillating pair-correlation function is what one anticipates when 
the system is in solid phase while a non-oscillating one signatures a liquid 
phase. 



djj(R,t) 
di 



(H — E T )ip(R, t) 
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1.5 Brief outline of the thesis 



The microscopic formalism developed in the first part of this thesis can 
accurately treat a weakly interacting Bose-Einstein condensate and a non- 
condensed thermal atoms for a range of temperatures below the 2D critical 
transition temperature T c . In the same spirit, I assume the semi-classical 
Two-fluid model to be valid in evaluating the equilibrium property of the 
condensate fraction and the density profile of bosonic gas in this limit. In 
this model the density profiles are calculated self-consistently by ensuring 
that the choice of effective interaction potential is introduced in a consistent 
manner. It is an extremely important point since the introduction of inaccu- 
rate effective potential costs heavily in the convergence of the numerical solu- 
tion. Another crucial fact to consider is the scattering property of a strictly 
2D Bose gas, where the s-wave scattering length has become larger than the 
axial thickness of the cloud changes dramatically. In Sec. 2.1 I will illus- 
trate how a semi-analytical expression for the many-body T-matrix elements 
corresponding to the condensate-condensate and to the condensate-thermal 
cloud coupling parameters (g% and gi, respectively) effects the equilibrium 
properties of the Two-fluid model at increasing temperature. 

Sec. 2.2 summarizes the Feshbach formalism framework with application 
to the 2D Bose gas. Here I examine the tuning of two-body collisions by 
means of a magnetic field in a strictly 2D condensate inside a harmonic trap 
at zero temperature. By modifying the coupling strength one can control 
the dimensionless physical parameter 77 which is the ratio between the mean- 
field interaction energy and the harmonic-oscillator level spacing |^T] . In 
the case of repulsive interactions (rj > 0), starting from a weakly interacting 
regime [r] <C 1) and up to a strongly interacting one (rj ^> 1), I evaluate the 
equilibrium density profile as a function of rj by solving numerically a non- 
linear Schrodinger equation by a split-step Crank-Nicholson discretization 
scheme. In the case of attractive interactions (rj < 0) one can give an upper 
bound for the number of bosons which can condense without the condensate 
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collapsing. At variance from the 3D gas (see for instance [JOl EJ), I have 
found that in 2D this critical number does not depend explicitly on the 
strength of the harmonic confinement. 

In Sec. 3.1, I consider a 2D rotating condensate at finite temperature and 
evaluate the particle density profiles and the energy of a vortex within a 
strictly 2D model for the boson-boson coupling. This is appropriate to a 
situation in which the s-wave scattering length starts to exceed the vertical 
confinement length. Previous work has established that the dimensionality 
of the scattering collisions strongly affects the equilibrium density profiles 
[751 IT§] and the process of free expansion of a BEC containing a vortex [HU] • 
A similar study of the density profiles of a rotating BEC in 3D geometry 
at finite temperature has been carried out by Mizushima et al. [HI], who 
also determined the location of various dynamical instabilities within the 
Bogoliubov-Popov theory. 

In Sec. 3.2, I will report the result of my investigation on the ground state 
properties of N v vortices treated as N Yukawa bosons trapped in the two- 
dimensional harmonic trap. A density functional theory (DFT) approach is 
developed to calculate the density profiles and the ground state energy of 
the vortex fluid within the Bogoliubov-de Gennes reference system. When 
a non-interacting auxiliary system is inferred in this model, we obtained a 
Kohn-Sham (KS) type equation. As a preliminary step, the density profile 
of the vortex fluid is calculated by solving this equation numerically. The 
method relies on the VMC data of Magro and Ceperley [71] on the homoge- 
neous ground state energies. The work in this thesis also pave way to many 
unresolved open questions pertaining to the physics of freezing and that of 
the strongly correlated incompressible liquid which has an overlap with the 
Fractional Quantum Hall Effect (FQHE). 
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Scattering properties for 2D 
Bose gas 

2.1 Coupling parameters in a mixture of con- 
densed and thermal-bosons 

The scattering property of bosons in a 2D flat geometry changes dramati- 
cally compared to the 3D one. The T-matrix for two-body collisions in vacuo 
at low momenta and energy, which should be used to obtain the collisional 
coupling parameter to the lowest order in the particle density, vanishes in the 
strictly 2D limit ^21 HH] • It is then necessary to evaluate the scattering pro- 
cesses between pairs of Bose particles by taking into account the presence of 
a condensate and a thermal cloud through a many-body T-matrix formalism 

inaisniiHi. 

This formalism has already been used in a number of studies of low di- 
mensional Bose gases, dealing in particular with phase fluctuations and the 
Kosterlitz-Thouless transition in a variational approach |52*1 183| l3*?l I85j . with 
a mean-field evaluation of the breathing-mode frequency in a trapped ID gas 
|8fij . and with the equilibrium density profile of a 2D condensate (HUES]- I n 
the limit of zero temperature the condensate-condensate coupling parameter 
has been related to the two-body T-matrix by considering that, when two 
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particles in the condensate collide at zero momentum, they both require an 
energy equal to the chemical potential /i to be excited out of the conden- 
sate [HI]. With increasing temperature the population of the excited states 
becomes non-negligible and a microscopic theory which also takes into ac- 
count the depletion of the condensate is required |87j . An additional coupling 
parameter to describe the scattering processes between an atom in the con- 
densate and an atom in the thermal cloud has been introduced by Stoof and 
co-workers [37| through the two-body T-matrix at energy — \i. 

The scattering processes between pairs of atoms in a gas consisting of 
a Bose-Einstein condensate cloud and a thermal cloud are described by 
the many-body T-matrix T MB (E) as a function of the energy E. The 
coupling parameters are given by the matrix elements (k.'\T MB (E) |k) = 
T MS (k, k', K; E), taken in the limit of zero energy and momenta. Here k 
and k' are the incoming and outgoing relative momenta of the pair of center- 
of-mass momentum K. I will consider only the condensate-condensate and 
condensate-thermal cloud couplings and neglect the scattering between ther- 
mally excited atoms in the present case of a dilute gas. 

Before discussing the many-body T-matrix, however, we shall first recall 
the behaviour of the two-body T-matrix that describes collisions between 
pairs of particles in vacuo. We shall then discuss how the two T-matrices are 
related in the appropriate limit for a two-fluid system. 



2.1.1 The two-body T-matrix 

The two-body T-matrix is the solution of the Lippmann-Schwinger equation, 

(k'\T 2B (E)\k) = (k'|y(|r 1 -r 2 |)|k) + ^(k'|y(|r 1 -r 2 |)|q) 

q 

x^^(q|T 2B (E)|k) , (2.1) 

with V(|ri — r2|) being the interparticle potential. The center-of-mass energy 
for the pair of atoms is E and each atom of mass m has a single-particle 
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excitation energy e q = h 2 q 2 /2m, as the collision takes place in free space. 
In Eq. (J2.1j) the collisions are described by a single-loop interaction between 
the two atoms plus contributions involving all possible transition routes from 
state |k) to state |k') iw intermediate states |q). 

In the case of the hard-disk potential of strength Vq = A-irh 2 /m and in the 
dilute limit (ka, k'a <1), the solution of Eq. ()2.1j) is [37| 

„,,„,,„.-.„ . 4iih 2 /m 

k' T (E) k « 2.2 

In ym l j Ema A \ 

In this case the T-matrix is independent of the momenta, but depends on the 
logarithm of the energy and vanishes as E approaches zero. Therefore, the 
presence of the surrounding gas must be taken into account in the collisional 
processes. Following the proposal of Morgan et al. [HEj, this is done by 
setting E = —2/j, for the mutual scattering of two condensate particles, this 
being the energy required for them to be excited out of the condensate. By 
a similar argument Al Khawaja et al. [HZ| set E = —\i for the scattering 
between a boson in the condensate and boson in the thermal cloud. 



2.1.2 The many-body T-matrix 

According to the above argument, a first approximation for the many-body 
T-matrix at zero momenta and energy is 

Tf>, 0, 0; 0) = T 2 >, 0, 0; -n/i) (2.3) 

with n = 1 for condensate-thermal cloud scattering and n = 2 for condensate- 
condensate scattering. Further many-body effects can enter the scattering 
processes attended by the presence of a condensate from attributing a Bo- 
goliubov spectrum to the intermediate states [B2]- At low momenta the 
many-body T-matrix is the solution of the integral equation 

(k'| T MB (E) |k) = (k'|V(|r 1 -r 2 |)|k) + ^(k'|V(|r 1 -r 2 |)|q) 

q 

31 



Coupling parameters in a mixture of condensed and thermal-bosons 



where Huj c 



\i are the Bogoliubov excitation energies in the Hartree 



approximation \i <C e q and iV(u; q ) ~ exp[/3(e q + /i)] — l" 1 are the corre- 
sponding population factors. For a contact interaction potential of Eq. ()2.4|) 
yields 

J_ _ 1+ 2JV(w q ) 



T M "(0,0,0;E) 



E - 2frw r 



(2.5) 



As discussed by Stoof and Bijlsma [S21ES1- Eq- ()2.5|) is still affected by an 
ultraviolet divergence, which can be remedied by replacing Vo in favor of the 
two-body T-matrix. The final result is 

iVK) 



T n MB (0,0,0;0)=T^(0,0,0;-n/.) 



2B, 



1+T 2B (0,0,0;-^)^ 



hhJ r 



(2.6) 



with T 2B (0, 0, 0; -n/i) given by Eq.fl22J with ^ 



-n/i. 



2.1.3 Calculation of coupling parameters 

The sum in the RHS of Eq. ([2.6j) can be evaluated analytically by replacing 
the sum over intermediate states by an integral over momentum. Setting 



N(u q 

hjJr, 



m 
2^ 



J2Ei(-sPn) 



(2.7) 



q ^ s=l 

where Ei(x) is the exponent-integral function. In the asymptotic low-temperature 
regime 3> 1) this thermal-population term can be approximated by 



J2Ei(-sPti -> ln[l -exp(-/^)]. 



(2. 



s=l 



A numerical illustration is given in Fig. (|2.1|) for values of the system pa- 
rameters appropriate to the experiments on 23 Na by Gorlitz et al. [H| (see 
Sec. (JUl). 



In summary, the coupling parameters in the 2D Bose gas are given by 

Auk 2 /m 



MB 
•J II 



In \4h 2 /(nmfia 2 )\ - 2 J2T=i Ei(-spn) ' 
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Figure 2.1: The correction term A(j3 ^)=J2^Li Ei{— s (5 ji) from excited-state 
occupancy (dashed line) and its approximate form Eq. ()2.8)) (solid line) as 
functions of j3fi. In the inset a zoom of the region 0.2 < /3fi < 1.2 is shown. 



with n=2 for collision between pairs of condensate atoms and n=l for colli- 
sions between an atom in the condensate and thermally excited atom. These 
parameters depend on temperature both through the chemical potential and 
through the excited-state population factor given in Eq. (|2.7j) and asymp- 
totically approximated by Eq. (j2.8|) . If this population factor is dropped, one 
obtains from Eq. (j2.3)l the corresponding "two-body" coupling parameters as 



Anh 2 /m 



In \4h 2 / (nm^ia 2 



(2.10) 



In the next Section I evaluate the chemical potential and hence the coupling 
parameters in Eqs. ()2.9j) and ()2.10j) through a self-consistent evaluation of 
the density profiles in a two fluid-model. 
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2.1 A Equilibrium properties in the two-fluid model 

Let n c (r) and n nc (r) be the particle density profiles for the condensate and 
for the thermal cloud in a 2D Bose gas which is radially confined inside a 

2 2 

isotropic planar trap described by the external potential V ext (r) = — j — . 
The two-fluid model [SHI EHj combines a solution of the Gross-Pitaevskii 
equation for the condensate with a Hartree-Fock model of the thermal cloud, 
which is treated as an ideal gas subject to an effective potential V e ff(r). In 
the present case, 



As already noted, we are neglecting the collisions between pairs of bosons 
belonging to the thermal cloud. 

It is well known that a full numerical solution of the Gross-Pitaevskii equa- 
tion Eq. (jl.28j) can be avoided when the kinetic energy term in it can be 
neglected |H)]. I postpone the evaluation of a full scale numerical solution 
of the GPE equation for the coming Sec. ()2.2|) on Feshbach resonances and 
in chapter 3, where the role of the kinetic energy become crucial in vortex 
calculation. Since our goal here is to analyse the effect of the interactions 
on the equilibrium properties of the system, we use the Thomas-Fermi ap- 
proximation in which kinetic energy is neglected. This yields a condensate 
density of the form 



n c (r) = (l/<? 2 )[/i - V ext (r) - 2g 1 n nc (r)]9( l i - V ext (r) - 2 9l n nc (r)). (2.12) 



where 9{x)= for x < otherwise 8(x)= 1 if x > 0. Since I am interested 
in examining qualitative behaviours rather than in attaining high numerical 
accuracy, I have adopted Eq. (j2.12| ) for the condensate density and found that 
discontinuities occurring in the density profiles at the Thomas-Fermi radius 
can be eliminated by the simple expedient of introducing a momentum cut- 
off in the expression for n nc {r) of Eq. ()1.30|) . That is, we calculate n nc {r) 
from 



Veff{r) = V ext {r) + 2 gi n, 



(2.11) 



rn 



Po 



,2 1 




2irh 2 f3 



In 1 - exp[P(fi - V e ff(r) 



m 



)] 



(2.13) 
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where we take po = ^Jvng\n nc (HD]. The model is then evaluated by solving 
self-consistently Eqs. (I2.11j) - (j2.13)) together with the condition that the areal 
integral of n c (r)+n nc (r) should be equal to the total number N of particles. 

The validity of the present model is as follows. A mean-field treatment is 
valid when the diluteness condition n c a 2 -C 1 holds and if the temperature 
of the gas is outside the critical region. With regard to the thermal cloud, 
no significant differences have been found between the predictions of the 
Hartree-Fock and Popov approximations in the regime n c a 2 <C 1, except at 
very low temperature where the thermal cloud is becoming negligible P?T] . 

2.1.5 Coupling strength and its effect on equilibrium 
properties 

For a numerical illustration, we have taken the values of particle number, the 
radial trap frequency, and the scattering length as appropriate for 23 Na atoms 
in the experiment of Gorlitz et n/.|T7j (N = 5 x 10 5 , u± = 188.4 Hz, a = 2.8 
nm). Whereas in their experiment collisions are in 3D regime, we focus on a 
strictly 2D regime that could be reached experimentally by increasing either 
the trap anisotropy parameter or the scattering length. 

I have compared the temperature dependence of the coupling parameters 
g 2B (long dashes for n=2 and short dashes for n=l) with that of the coupling 
parameters g^f B (full and dotted lines, respectively) in Fig. (|2.2[) . It is evi- 
dent that the many-body screening of the interactions due to the occupancy 
of excited states is quite large and rapidly increasing with temperature. 

Such many-body screening has, however, very little effect on equillibrium 
properties of the gas for our choice of system parameters. Figure ()2.3j) reports 
the condensate fraction N /N for the g^f model as a function of temperature, 
in comparison with that of an ideal Bose gas at the same values of the system 
parameters. As is well known, the transition temperature and the condensate 
fraction are lowered by the interactions. However, the g 2B model gives results 
that are practically the same as the g^f B one. 
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Figure 2.2: Interaction strengths (in units of h 2 /m) as functions of temper- 
ature T (in units of ideal-gas critical temperature To). 
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Figure 2.3: Condensate fraction Nq/N as a function of temperature T (in 
units of To) from the g^ B model and (dashed line) compared with the non- 
interacting gas (solid line). 
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The evolution of the density profiles for the condensate and for the thermal 
cloud with increasing temperature from near absolute zero to the critical 
temperature T of the ideal gas (T = (y/E/ir)fiw\/N jH21) are depicted in 
Figure ()2.4|) . The results in Figure ()2.4b ) are in good agreement with those 
of Tanatar et al. [78] . except that the tails of that profile are missed in 
the Thomas-Fermi approximation. Again the screening of collision from the 
occupancy of the excited states is very small and becomes barely visible at 
T « 0.75 T . 

We can deduce that the use of Eq. ()2.3|) to describe the many-body effects 
in two-body scattering processes in 2D Bose-condensed gas appears to be 
very good in regard to equilibrium properties at large values of the particle 
number. A decrease in the number of particles lowers the chemical potential 
and may lead to observable effects for N f=s 10 3 . 
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Figure 2.4: Density profiles of the condensate (solid line) and the ther- 
mal cloud (dashed line) in the g^ IB model (in units of N/a\ with = 
(h/muj) 1 / 2 ) as functions of radial distance r (in units of a/j Q ) at various 
values of the temperature ( T/T =0,0.3,0.5,0.8 and 0.95 from (a) to (e)). 
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2.2 Feshbach resonance in a 2D Bose gas 

The 3D s-wave scattering length can be tuned by means of external electric 
or magnetic fields HUH US] • Within this process a resonance can appear 
at a particular value of the magnetic field and allows control of the coupling 
strength ranging from positive (repulsive) to negative (attractive) values. 
One may study new dynamical effects of collective oscillations in large con- 
densates jHH] and the dynamics of collapsing and exploding condensates [§Jj] . 
Calculations of such Feshbach resonances have been given for 3D gases of 7 Li 
[HI!, 23 Na (Eg , 39k and 4lK US], 85 Rb and 87 Rb [EJ, and 133 Cs jHEJ. Experi- 
mental observations have been reported for 87 Rb 85 Rb HTM HTHl ITTT2] . 
23 Na |lU3j . and 133 Cs |lU4j . Feshbach resonances in quasi-2D atomic gases 
have been discussed by Wouters et al. |lU5j . who have predicted a shift in 
the position of the resonance from the squeezing of the confinement. 

2.2.1 Overview on scattering theory 

We start by recalling some well-known facts about a collision between two 
particles with internal degrees of freedom governed by an internal Hamilto- 



where |a) and e a are the single-atom eigenstate and its corresponding eigen- 
value. For example the internal Hamiltonian can describe the hyperfine in- 
teraction between nuclear and electronic spins and their Zeeman coupling 
to an external magnetic field B. A two-body collision is described by the 
Hamiltonian 



in the center-of-mass system, with p the relative momentum and m r the 
reduced mass (see for instance ^06 ). The Hamiltonian in Eq. ()2.15j) is 
the sum of a part with eigenstates \a(3) ("channels") tending asymptotically 
to a symmetrized product of separate-atom internal states |a) and and 



nian H 



H int \a) 




(2.14) 




(2.15) 
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of a finite-range interaction V(r) which couples the channels. The energy 
associated with an eigenstate \a/3) is E a) p = E a +e f 3+ti 2 kl l ^/2m r , where hk^p 
is the relative momentum of the two incoming particles. The asymptotic 
scattering state ^^(r) can be expanded on all scattering channels in the 
form 

tt a/J (r) = exp(zk a/3 • r) \a/3) + ^ f*f (Kp, k^) ^^' ' ^ W P') , 

(2.16) 

where the outgoing channel is denoted by \ct' (3') and the scattering am- 
plitude f^f is proportional to the T-matrix element (k' Q//3 , | T |k Qj a). The 
asymptotic magnitude of the momentum for the \a'/3') channel is k a ip> = 
a/3 — e a' — £p'), from conservation of total momentum and energy 
in the scattering process. The channel is said to be open (closed) if the 
momentum takes positive real (imaginary) value. 

The whole Hilbert space describing the spatial and spin degrees of free- 
dom is now divided into a subspace P for open channels and a subspace Q 
including all closed channels |lU7j . The state vector and the Schodinger 
equation H = E |\&) are projected onto the two subspaces by means of 
projection operators P and Q satisfying the conditions P + Q = 1 and PQ = 
0. The formal solution of the coupled equations obeyed by the two projected 
components is 

\9 Q ) = {E- Hqq + iSy'HQp \9 P ) (2.17) 

and 

(E-Hpp-H' PP ) |W F ) = 0, (2.18) 
where H PP = PHP etc., 

H'pp = H PQ (E - H QQ + i5)- 1 Hq P , (2.19) 

and 5 is a positive infinitesimal. The term H' PP in Eq. ()2.18|) describes the 
Feshbach resonances. It represents an effective non-local, retarded interaction 
in the P subspace due to transitions to the Q subspace and back. The 
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Hamiltonian Hpp + H' PP gives the effective atom- atom interaction in the 
open-channel subspace, the effects due to the existence of bound states being 
taken into account through the term H' pp . 



2.2.2 Feshbach resonances in the 2D coupling 

A Feshbach resonance results when true bound states belonging to the closed- 
channel subspace match the energy of open channels: transient transitions 
may then be possible during the collision process [SH]. Setting H PP = H + Ui 
where Hq includes the kinetic energy of relative motion and the internal 
Hamiltonian, the scattering matrix elements in the P subspace between 
plane-wave states with relative momenta k and k' read 

(k'| T |k) ~ (k'| T x |k) + (Oi,+k'| H' PP |fii, + k) • (2-20) 

Here 71 = U\ + \J\G§Tx represents the scattering process if Q subspace is 
neglected and = (1 — GqUi)~ 1 is the wave operator due to the interaction 
potential U 1} with G = (E — H G )~ l [gj. In Eq. (j2~20jl the effect of closed 
channels has been included at first order in H' PP . In the limit of zero relative 
velocity Eq. ()2.20|) leads to an expression for the coupling strength, 



9 = 9 



| (lpn\H Q p\jjo) | 2 ^ 

E t h — E<- 



where the sum is over all the states \ip n ) in Q subspace and we have set 
|Jli i+ k) ~ |fii j+ k') = |^ )- Here, g is the coupling strength estimated in the 
absence of closed channels and E t h is the threshold energy for the state |^o) 
at vanishing kinetic energy. Finally, when the threshold energy is close to 
the energy E res of one specific bound state \ip res ), this contribution will be 
dominant and the contributions from all other states may be included in g 
through an effective non-resonant scattering length a nr . For a strictly 2D 
condensate we get 

4irH 2 /m + | {^ res \H QP \iJo} | 2 ^ ^ 



In \Ahj \ima? nr \ E th — E, 
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where /i is the chemical potential of the condensate. Equation (|2.22J) is our 
main result. The non-resonant first term has been calculated from the two- 
body T-matrix for a binary collision occurring at low momenta and energy 
in the presence of a condensate, which fixes the energy of each colliding atom 
at the chemical potential [EEllISj- The resonant second term can be tuned by 
exploiting the dependence of the energy denominator on external parameters 
and in particular on the strength of the applied magnetic field. Explicitly, if 
the energy denominator vanishes at a value B of the field, then by expansion 
around this value we have 

E th - E res « (m res -m a - m /3 )(B - B ) (2.23) 

where m a ,p = —de a ^/dB are the magnetic moments of the two colliding 
atoms in the open channel and m res = —8E res /dB is the magnetic moment 
of the molecular bound state. Equation ()2.22|) then yields 

- 4rrh2 ( 1 AB \ 

9 ~ ~rn~ \\n\4h/fima 2 nr \ + B — B ) ' ^ > 

where the width parameter AB is defined as 

AB = J±_\^res\H QP \^)\\ 

Anh 2 m res — m a — 

Higher-order terms in H' PP will lead to a broadening of the resonance. How- 
ever, the width of resonant states close to the threshold energy in the open 
channel is usually small because of the low density of states, and Feshbach 
resonances are usually very sharp jH]. 

In the next Section I will analyze the consequences of tuning the coupling 
strength on the equilibrium state of a 2D condensate in a harmonic trap. 



2.2.3 Equilibrium properties as functions of the cou- 
pling 

We consider a Bose gas with a fixed number N of particles in a 2D harmonic 
potential V ext (r) = muj 2 r 2 /2. At zero temperature the condensate wave 
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function \l/(r) is described by the mean-field Gross-Pitaevskii equation 



h 2 1 d 



2m r dr \ dr 



d 

) + V ext (r) + gn c (r) 



v[/(r) = /i\l/(r) 



(2.26) 



where n c (r) = |\l/(r)| 2 is the density of the condensate. We introduce the 
dimensionless coordinate x = r/a^o, where a^ = (h/mto) 1 ^ 2 is the harmonic 
oscillator length, and the dimensionless chemical potential /i' = 2fi/hw. Thus 
Eq. ()2.2(i|) can be written in rescaled form as 

2 



d 2 



1 al 



dx 2 x dx 



^>(x) + x 2 ^(x) + 2?7 



*(0) 



= ii'V(x) , (2.27) 



where rj — g \^f(0)\ 2 /huj is a dimensionless coupling parameter given by the 
ratio between the mean-field interaction energy per particle and the oscillator 
level spacing. The normalization condition on the wave-function is 



N = a 



ho 



xdx\^> (x) 



(2.28) 



Equations (|2.27|) and (J2.28J) need to be solved self-consistently. 



2.2.4 Repulsive coupling (77 > 0) 

In the case of repulsive interactions Eq. (|2.27|) always admits a solution: 
the condensate exists and is stable in the 2D harmonic trap. Of course, far 
away from the resonance the dimensionless coupling parameter rj is primar- 
ily determined by the non-resonant scattering length. On approaching the 
resonance by increasing B towards B , the coupling parameter becomes very 
small for B ~ B$ — AB/(ln \&K/ /xma^ r |): the non-linear term in Eq. ()2.27|) 
becomes small and the condensate wave function is close to a Gaussian pro- 
file. On further increasing B the gas enters the strong-coupling regime and 
in the limit 77 ^> 1 the kinetic energy term in Eq. ()2.27|) and the discrete 
structure of the trap levels become irrelevant. The profile then approaches 
a Thomas- Fermi form |3*3*| l4*H] . The shape of the condensate density pro- 
file, obtained by the numerical solution of Eq. ()2.27|) at several values of 
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Figure 2.5: Scaled density profile |\l/(x)| 2 /|\l/(0)| 2 as a function of the dimen- 
sionless radial distance x = r/aho for various values of the coupling parameter 
t]. Dashed curves are Thomas-Fermi profiles and the dotted curve is the free- 
gas profile. 

the coupling parameter rj, is depicted in Fig. 12.51 The Figure illustrates 
comparisons with the free-gas and Thomas-Fermi profiles, according to the 
discussion given just above. 

2.2.5 Attractive coupling (77 < 0) 

On going through the resonance the Thomas-Fermi-like density profile will 
suddenly shrink and collapse. Indeed, in the strong-coupling attractive regime 
the kinetic energy does not suffice to compensate for the negative mean-field 
energy term and Eq. (|2.26J) does not admit physical solutions. A condensate 
with a limited number of particles can only exist in the weakly attractive 
regime. An upper bound for the critical number N c of bosons in a conden- 
sate with weakly attractive interactions can be estimated by comparing the 
kinetic energy and the mean-field energy at the centre of the trap. This leads 
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to the inequality 

1 + 2rj > (2.29) 

for condensate stability, where in the calculation of the kinetic energy we 
have used the free-gas wave function. Using the expression of rj we find from 
Eq. (|2~29j) 



1 AB 



-i 



\n\4h/fima 2 nr \ B — B ' 
At variance from the trapped 3D gas jlH], the critical number N c does not 
depend on the harmonic oscillator length, the dependence being cancelled 
by the fact that in 2D the mean-field energy scales as ~ like the kinetic 
energy. Let us remark that the condensate in the case of attractive coupling 
can be metastable even if the number of condensed bosons is lower than N c . 
A three-body loss term should be included in Eq. ()2.2fij) in order to calculate 
the decay time and to describe the time evolution of such a condensate, as 
was done for the 3D Bose gas [J3J 021 EH] and very recently for a 3D boson- 
fermion mixture |lU8j . An evaluation of three-body losses in the 2D Bose 
gas and of the time evolution of a 2D condensate in the weakly attractive 
coupling regime may be considered in the future. 
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Chapter 3 

Rotating 2D trapped Bose gas 



3.1 Single vortex in a rotating 2D Bose gas 



Understanding the behaviour of vortices in the 2D regime is important, since 
they reflect the superfluid nature of the condensate |1U91 II 101 l§T] and are 
expected to play a role in the transition from the superfluid to the normal 
state ^Hj- An experimental method for the creation of quantized vortices 
in a trapped BEC has made use of an "optical spoon" j3U llllj . whereby 
the condensate in an elongated trap is set into rotation by stirring with 
laser beams. This experiment on a confined boson gas is conceptually the 
analogue of the rotating bucket experiment on bulk superfluid Helium, with 
the difference that the thermal cloud is inhomogeneously distributed and lies 
mostly outside the condensate cloud. A quantized vortex first appears in the 
condensate at a critical angular frequency of stirring which corresponds to an 
instability of the vortex-free state. A lower bound for the critical frequency 
can be assessed from the energy of a vortex, defined as the difference in 
the internal energy of the gas with and without a vortex. Observation of a 
vortex in a trapped gas is difficult since the vortex core is small in comparison 
to the size of the boson cloud, but the size of the core increases during free 
expansion and indeed vortices were first observed experimentally by releasing 
the trap and allowing ballistic expansion of the cloud [q"4*] . 
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3.1.1 GP model that accomodates a single vortex 

The BEC is subject to an anisotropic harmonic confinement characterized 
by the radial trap frequency u± and by the axial frequency Xu± with A ^> 
1. Motions along the z direction are suppressed and the condensate wave 
function is determined by a 2D equation of motion in the {x, y} plane. 

The order parameter for a 2D condensate accommodating a quantized 
vortex state of angular momentum Uk per particle is written as ^(r) = 
i/j(r) exp(w0), with <fi the azimuthal angle. The wave function ^f(r) then 
obeys the nonlinear Schrodinger equation (NLSE) 



H 2 2 . fr 2 K 2 1 



*(r) = ii$(r) (3.1) 



similar to Eq. (jl.28j) but in addition it includes the contribution 2gin nc (r) 
of the thermal clouds. The 2D coupling parameters gj, with j = 2 for 
condensate-condensate repulsions and j — 1 for condensate-noncondensate 
repulsions, are given by 



Ann 2 Ira 
yj ln|4/i 2 /(jm/ia 2NI ' 



In Eq. ()3.2|) above we have omitted a term due to thermal excitations [75] . 
which is negligible in the temperature range of present interest (T < 0.5T C , 
with T c the critical temperature). 

In this model the atoms in the thermal cloud are not put directly into 
rotation, but feel the rotating condensate through the mean-field interactions. 
In the Hartree-Fock approximation j2H] the thermal cloud is treated as an 
ideal gas subject to the effective potential V e ff(r) written as in Eq. (|2.11j) . 
The density distribution of the thermal cloud is then given by Eq. (j2.1Hj) 
in which the momentum cut-off po in Eq. (|2.1H|) is equivalent to adding the 
term 2gin nc to the effective potential in Eq. (j2.11|) . 

We solve self-consistently the coupled Eqs. (13. lj) and (j2.13j) together with 
the condition that the areal integral of n c (r)+n nc (r) is equal to the total 



48 



Single vortex in a rotating 2D Bose gas 



number N of particles. The differential equation ()3.fj) is solved iteratively 
by discretization, using a two-step Crank-Nicholson scheme |112j . In Sec. 
(J3. 1.3)1 I will also compare the results with those obtained in the Thomas- 
Fermi approximation by dropping the radial kinetic energy term in Eq. 1)3.1)1 . 



3.1.2 The energy of a vortex 

The critical angular velocity measured in experiments where vortex nucle- 
ation occurs from a dynamical instability |3] depends strongly on the shape 
of the perturbation. However, a lower bound for the angular velocity required 
to produce a single-vortex state can be estimated once the energies of the 
states with and without the vortex are known. Since the angular momentum 
per particle is Hk, the critical (thermodynamic) angular velocity is given by 

«. = ^|=*. (3.3) 

This expression follows by equating the energy of the vortex state in the 
rotating frame, that is E K — tt c L z , to the energy E k= q of the vortex- free 
state. 

In the noninteracting case at zero temperature, the energy difference per 
particle is simply Iiku)±, so that fl c is just the trap frequency in the {x,y} 
plane. For the interacting gas at zero temperature in the Thomas-Fermi 
approximation, Eq. ()3.3)1 reduces to the expression 1109) 

nr(0) ^ ln (^). 

Here R = {2n/rmj\) 1 / 2 and £ = Rfrw ± /2fi are the Thomas- Fermi radius and 
the healing length respectively. 

In the general case of an interacting gas at finite temperature, we have to 
evaluate numerically the total energy as the sum of four terms 



E K — -Ekin,c + ^trap + ^int + -Ekin,T • (3.5) 
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These terms are the kinetic energy of the condensate 

£^ = /^<r)(-£v> + |£l)*(r), (3.6) 
the energy of confinement 

-Strap = \ rmj '±- J d 2 rr 2 [n c (r) +n nc (r)] , (3.7) 
the interaction energy 

E int = ^J d 2 r[g 2 n 2 c (r) + A 9l n c {r)n nc {r) + 2 9l n 2 nc {r)] , (3.8) 
and the kinetic energy of the thermal cloud 

= J d 2 r J {exp[(3(p 2 /2m + V eff (r) - /,)] - l}" 1 . (3.9) 

In Sec. f!3.1.3|) I compare the angular frequency obtained by the calculation 
of the total energy of the gas with and without a vortex with those obtained 
from a Thomas-Fermi calculation and from extrapolating Eq. (J3.4)) at finite 
temperature through the temperature dependence of the chemical potential. 

3.1.3 Numerical result for single vortex 

For a numerical illustration we have taken k = 1 and chosen values values 
of particle number, the radial trap frequency, and the scattering length as 
appropriate for 23 Na atoms in the experiment of Gorlitz et a/.^Tj (N = 
5 x 10 5 , u>± = 188.4Hz, a = 2.8 nm). I implicitly assume, however, that the 
trap has been axially squeezed to reach the strictly 2D scattering regime. 
The density profiles obtained from the solution of Eqs. (J3.1)) for the gas at 
three different values of the temperature (in units of the critical temperature 
T c = (\fON liik^jhhj^ of the ideal Bose gas) are shown in Figs. 13. II and 13.21 
In the absence of a vortex (Fig. 13. lj) . the main point to notice is that the 
growth of the thermal cloud exerts an increasing repulsion on the outer parts 
of the condensate, constricting it towards the central region of the trap. 
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Figure 3.1: Density profiles n(r) for the condensate and the thermal cloud 
(in units of a 2 ho /N ', with a^o = \fh]muJ±) versus radial distance r (in units 
of aho) at temperature T/T c = 0.05 (full line), 0.25 (long-dashed line), and 
0.50 (short-dashed line). 



This effect is seen only when the condensate is treated by the NLSE and 
persists in the presence of a vortex (Fig. 13. 2 J) , but as we shall see below, is 
missed in the Thomas-Fermi approximation where the radial kinetic energy 
of the condensate is neglected. In addition, the thermal cloud penetrates 
the core of the vortex and enhances the expulsion of the condensate from 
the core region (inset of Fig. 13. 2|) . in a manner which again is governed 
in its details by the radial kinetic energy term in the NLSE. The profiles 
obtained for a BEC containing a vortex at T = 0.5 T c by the full numerical 
calculation using the NLSE are compared in Fig. 13.31 with those obtained 
in the Thomas-Fermi approximation. The constriction of the condensate in 
its outer parts and its expulsion from the core region by the thermal cloud 
are clearly underestimated in the Thomas-Fermi theory. Finally, Fig. 13.41 
reports our results for the energy of the vortex as a function of T/T c . The 
inaccuracies arising in the density profiles from the Thomas- Fermi treatment 
of the interplay between the condensate and the thermal cloud clearly lead 
to large errors in the estimation of Q C (T). On the other hand the simple 
expression given in Eq. (|3.4|) . with the two alternatives of using in it the 
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Figure 3.2: Density profiles for the condensate and for the thermal cloud in 
the presence of a vortex. Units and symbols are as in Fig. 13.11 The inset 
shows an enlarged view of the profiles near the center of the trap. 
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Figure 3.3: Density profiles for a BEC containing a vortex at T/T c = 0.50 
in the full calculation using the NLSE (full lines) and in the Thomas-Fermi 
approximation (dashed lines). The units are as in Fig. 13.11 The inset shows 
an enlarged view of the profiles near the center of the trap. 
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Figure 3.4: The energy of a vortex, expressed as the thermodynamic critical 
frequency fl c in units of the radial trap frequency u>±, as a function of T/T c . 
The results from the full calculation using the NLSE (full line) are compared 
with those obtained from Eq. ()3.4j) with the corresponding values of /x(T) 
(dotted line) and with the Thomas- Fermi values of /i(T) (short-dashed line). 
The long-dashed line shows the results obtained in a calculation using the 
Thomas-Fermi theory. 
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chemical potential fi(T) from the Thomas-Fermi approximation or from the 
full calculation using the NLSE, gives a reasonable account of the vortex 
energy up to T ~ 0.5 T c . 

3.2 Density Functional Theory for many- vortex 
problem 

3.2.1 Introduction 

Nelson et al. [73] have shown that the statistical mechanics of the flux- 
line lattice (FLL) of high-T c superconductors can be studied through an 
appropriate mapping onto the 2D Yukawa boson system interacting with 
potential Vy{r) = eK (r/a) where e is an energy scale and a is a length scale 
having the meaning of a screening length. The function Kq(x) being the 
modified Bessel function and the coefficient e attached to it is the coupling 
strength that scales the energy of interaction. 

The Hamiltonian of N bosons describing this system reads 

ff=-£A a V? + £#o(r y ). (3.10) 

i i<j 

Here i, j are particle indices and the de Boer dimensionless parameter is 
defined by A 2 = h 2 /2ma 2 e. Meanwhile the reduced density of the system is 
denoted by p = N/Aa 2 , where A is the surface area. These two dimensionless 
parameter characterize the system. 

On the other hand, the Hamiltonian for a system of N vortices mapped as 
a system of N boson interacting with Yukawa potential is written as 

F = _ s M2! v? + s _^_ A . o(rw/A) _ (311) 

% Kj 

The system described by the Hamiltonian ()3.11)) is equal to one described by 
Hamiltonian (l3~TU|) if we set a = X and A 2 = (27rA; B T) 2 )/2e^. 
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Figure 3.5: The phase diagram of Yukawa bosons reproduced from Magro 
and Ceperley [22! . The pluses are transition points computed with DMC 
The dashed lines at high densities is the scaling law: A ~ 0.04/y/p. 

Following this idea, Magro and Ceperley [221 EH] have performed a Dif- 
fusion Monte Carlo (DMC) and Variational Monte Carlo (VMC) numerical 
calculation to calibrate the ground state vortex properties for a homogeneous 
2D Yukawa boson system. Of particular interest to us is the phase diagram 
obtained for the data of (p, A) at transition points displayed in Fig ()H.5j) . 

They observed that for a particular region A > 0.09 the system is domi- 
nated by kinetic energy and does not crystallize. Below this threshold how- 
ever a peculiar behaviour of reentrant liquid has been observed. Which means 
that system is in the liquid phase at very low density as well at high density. 
The crystal melts on compression and expansion. Nevertheless this reentrant 
liquid behaviour has not been fully understood. 

Inspired by their work, I will investigate the ground state property of a sys- 
tem of N vortices mapped as N bosons interacting via 2D Yukawa potential 
V(x) = eK (x) confined in a 2D harmonic planar geometry |ma>^r 2 . As a 
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matter of fact, the Yukawa potential has an interesting combination of short 
range with soft core. For small x, Kq(x) diverges like — ln(x) while for large 
x it decays as exp(— x)j \fx. The logaritmic term represents the inter-vortex 
repulsion and the exponential decay signifies that vortices far apart do not 
feel each other. 



3.2.2 Density Functional Theory 

With the Yukawa potential we cannot use the similar approach used for the 
contact delta potential Vo<5(r — r') in deriving a dynamical equation of motion 
as described in chapter 1. Here we use the Density functional theory (DFT) 
which is originally based on the notion that for a many-electron system there 
is a one-to-one mapping between the external potential and the electron 
density: v ext (r) <-> p(r). In other words, the density is uniquely determined 
given a potential, and vice versa. All properties are therefore a functional of 
the density, because the density determines the potential, which determines 
the Hamiltonian, which determines the energy and the wave function . 

Following this train of reasoning, the inhomogeneous dilute system of N 
interacting bosons can be described within the second quantization language 

as 



H 



Ho + J dr^(r)V ext (r)^{r) 

\J dr f dr'ip\r)^(r')V(\r -r'\)iP(r')iP(r) 



H + V ext + V t 



hit 



(3.12) 



where H$ = j drip\r) — |^;V 2 — \x ip(r) and V(|r — r'|) is the inter-atomic 
interaction potential. The annihilation and creation field operators are de- 
noted by "0'(r) and ip(r') respectively and obey Bose-Einstein commutation 
relations: 



,^(r , )] = [V t (r),V t (r , )]=0. (3.13) 
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Let us denote the ground state of the system as \g) so the ground state 
energy is defined as E = (g\H \g) and the density as n(r) = (g \ ip^ip \g). The 
Hohenberg-Kohn (HK) theorem |113j guarantees that there exists a unique 
functional of the density, 

F[n(r)] = (g\H + V int \g) , (3.14) 

irrespective of the external potential. The theorem was originally proved for 
fermions but its generalization also covers bosons. Following HK, we can 
write the total energy functional of the system as following, 

E[n{r)} = F[n{r)} + J drV ext (r)n(r) (3.15) 

Determination of the ground state energy follows by imposing the stationary 
conditions 

SE\n(r)] 

W = (3 ' 16) 

In general, the ground state cannot be determined exactly, but one has to 
resort to the Kohn-Sham procedure |114j to introduce an accurate approxi- 
mation. The idea is to map the interacting system of interest to an auxiliary 
system. In order to describe the auxiliary system, let us first decompose the 
Bose quantum field operators ip(r) as a sum of spatially varying condensate 
$(r) and a fluctuation field operator if) : 

^(r) = $(r)+$(r) , ^(r) = $*(r) + ^ f (r) (3.17) 

The wave function $(r) is also known as the order parameter and is defined 
as the statistical average of the particle field operator $(r) =< >. On the 
other hand, the condensate density is defined in term of order parameter as 
n c (r) = |$(r)| 2 . In the well-known Bogoliubov approximation [23] for a dilute 
weakly interacting gas, in which almost all the atoms are Bose condensed, one 
keeps only the term up to quadratic in the non-condensate field operators 
ip and i\v . When the fluctuation field is comparably small to the spatial 
order parameter (in an average sense, <C one can comfortably utilise 
a Mean-field approach. Using the decomposition in Eq. (|3.17J) above and 
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applying the Bogoliubov approximation, one can replace the product of field 
operators (r)ip^ (r')ijj(r')ip(r) in the interaction energy term of Eq. (J3.12j) by 

^(r^VMr'Mr) = |$(r)| 2 |$(r')| 2 + 2|$(r)| 2 \<$>(r')ft(r') + $*(r')^(r') 

+ |$(r')| 2 ^ t (r)^(r) + $*(r')$(r)ft(r)ft(T f ) 
+ i$*(r / )$*(r)^(r)^(r / ) + ^(r')^>(r)^(r)^(r'j3.U 

So now we are in a position to define the Hamiltonian of the auxiliary system 
that includes the inter-particle potential defined above as 

h 2 



2m ^ ^ 



Hr) + V s + J dv^{v)V: xt {v)^{r). (3.19) 



where our definition of auxiliary interparticle potential reads 

Vs = I dr J dr'V{\r -r'\)\${r)\ 2 <f>(r')ft(r') + <$>*(r')4>{r') 

dr I dr'V{\r-r'\)[2^*(r')^(r)^(r)ij(r') 



+ <$>(r')®(r)ft(r)^(r') + $*(r')$*(r)^(r)^>(r')] 



(3.20) 



Thus for the ground state \g s ) we can define the unique ground state energy 
functional of the auxiliary system as 

h 2 



F s [n s (r)} = (g s \ J dr^(r) 



2m 



4(r)\9.) + <9.\Vs\g.) , (3.21) 



which does not depend on the auxiliary external potential V^ xt . Hence, the 
total energy functional of the auxiliary system can be written as 



E s [n s {r)] = F s [n s (r)} + / drn(r)^(r 



(3.22) 



in which E s [n s (r)] can be approximated to -E^r)] or in other word the 
density of the auxiliary system n s (r) is identical to the real system n(r) by 
choosing a proper choice of auxiliary external potential V^ xt . 

Since the non-interacting part of Eq. ()3.21)1 is equal to the non-interacting 
part of Eq. ()3.14)1 and comparing terms in Eq. (J3.2()j) with Eq. (J3.18j) we 
obtain the following functional relation: 

F[n(r)] = F s [n(r)] + J dr J dr V{r - r')n(r)n(r') + F xc [n(r)] . (3.23) 
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The second term Vh is called the Hartree-energy defined as 

V H = ^Jdr j dr'V{\r-r'\)n{r)n{r') (3.24) 

in which the total density n(r) includes a non-condensate local density de- 
fined as 

ra(r) =< ^(r)^(r) >= n(r) - |$(r)| 2 . (3.25) 

The last term in Eq. ()3.23j) represents the exchange-correlation energy 
E xc [n(r)] and includes all the contributions to the interaction energy be- 
yond mean field. Calculating the variational derivatives in Eq. (j3.16|) using 
Eq. (|3~23j) . one finds 

5F s \n(r)] , . 6F xc \n(r)} , , , 

where the Hartree field read 

V h (t) = J dr'V{\r-r'\)n(r'). (3.27) 

Performing a similar variational calculation on Eq. (J3.22|) we deduce that 
the density of the auxiliary system is identical to the actual system if 

V: xt {r) = VUO + V H (v) + . (3.28) 

3.2.3 Interacting auxiliary system : Bogoliubov ap- 
proach 

In general, if an interacting auxiliary system is chosen, one needs to do 
further analysis by using a Bogoliubov approach. We rewrite the auxiliary 
Hamiltonian Eq. (j3.19j) using the decomposition in Eq. (j3.17|) by 

H s = J dr$*(r)L$(r) + J dr^t(r)L$(r) + J dr$* (r)L^(r) 

+ j drip\r)L4)(r) + dr J dr'V Y (\Y - r'|)[2$*(r / )$(r)# t (r)^(r') 
+ $(r')$(r)$ t (r)$ t (r / ) + $*(r')$*(r)^(r)^(r')] (3.29) 
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where the operator L is defined by 

In order to diagonalize Eq. ([3.29)1 . we first eliminate the terms linear in ijj 
and i\r by requiring that $(r) satisfy the equation 

L$(r) = 0. (3.31) 

Hence our auxiliary Hamiltonian Eq. ()3.29j) reduces to 

H s = J drft(r)Li>(v) + j dr J dr'V{\r - r'|)$*(r / )$(r)^ t (r)^(r') 

+ \J dr J dr'V(\r-r'\)[$(r')§(r)ft(r)ft(r') 

+ $*(r')$*(r)^(r)^(r')] . (3.32) 

The quadratic expression given by this equation can be diagonalized by the 
usual Bogoliubov transformation [23] 

i 

^W^EKW^-^'W^']' ( 3 - 33 ) 

i 

where the quasiparticle operators a, and a\ satisfy Bose commutation rela- 
tions. The values of Ui and Vi are obtained by solving the following coupled 
Bogoliubov equations: 

Lui + j dr'T/(|r-r'|)$V)$( r ) w i( r/ ) 

+ Jdr'VQr- r'\)<f>(r')${r)vi{r') = E iUi (r) (3.34) 



and 



Lvi + J dr'V{\r-r'\)$*{r)${r')vi{r') 

+ J dr' V{\r - r'|)$*(r)$*(r)« i (r') = -E iVi {r) . (3.35) 
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Hence with the solution of U{ and Vi from the above coupled equations, the 
Hamiltonian in Eq. (J3.29|) can be shown to reduce to 

K ew = -J2 E * I dr\ Vj (r)\ 2 + J2EM*i, (3-36) 

i J i 

which describes a non-interacting gas of quasiparticles of energy Ej. The 
ground state expectation value of Eq. ()3.36|) is given by 

E s = -Y,E l I dv\ V] {v)\\ (3.37) 

i 

taking into consideration that the ground state \g s ) obeys a \g s ) = 0. Besides, 
using Eq. (J3.33J) . the local non-condensate density defined in Eq. ()3.25|) can 
be written as 

n(r') = (&|^(r)^(r) \g s ) = ^|^(r)| 2 . (3.38) 

i 

Thus performing a little algebra, the new auxiliary functional which includes 
depletion of the condensate is defined by 

F s [n(r)] = J dr$*(r) (-—- - $(r) ~J2 J dr l^( r )H^ + ^( r )] • 

(3.39) 

The key element in the above derivation is that one can take into account the 
depletion of the condensate as shown in Eq. ()3.38)1 . Using the the auxiliary 
system defined by the Hamiltonian Eq. (J3.19)) and the relation in Eq ()3.28j) . 
both the condensate density n c = |$(r)| 2 and also the total density n(r) can 
be computed. 

The difficulties in finding the solution for Eq. (|3.31|) . Eqs ()3.34j) and 
(J3.35J) lies in the fact that the exact functional dependence of the exchange- 
correlation potential w xc [n(r)] = 8F xc [n(r)]/5n(r) is unknown for most of the 
system of interest and thus one has to resort to approximations such as the 
Local Density Approximation (LDA). 
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3.2.4 Application of DFT to N vortex problem 

We are now in a position to analyse our problem of N vortices treated as N 
Yukawa bosons via the one-to-one mapping of Nelson and Seung [ZHj with the 
DFT tool summarized in the previous section within our approriate scaling. 
By choosing an auxiliary system of N non-interacting bosons (referring to 
Eq. (j3.21j) and Eq. (|3.28jl ). Eq. (j3.31j) then can be written more explicitly 
as satisfying the following equation: 



This equation is a sort of generalization of Gross-Pitaevskii equation for a 
dilute inhomogeneous gas at T = 0. Since we have chosen a non-interacting 
auxiliary system, we expect that at T = all the atoms are Bose condensed 
as their total density n(r) is approximated by condensate density n c (r) = 
|$(r)| 2 . Setting F xc [n(r)} = one recovers a closed non-linear Schrodinger 
equation (NLSE). 

The available numerical data of Magro and Ceperley [13] and Strepparola 
et. al |115j permit one to obtain information on the homogeneous excess 
free energy /e^b 7 -] rather than the homogeneous exchange correlation energy 
f xc [n] . Thus it is much convenient to work with excess free energy defined 
as [TT6] 



In general, the excess correlation functional energy F ex [n(r)] in a DFT cal- 
culation is not known exactly. One can resort to approximations such as the 
Local density approximation (LDA) which reads, 



where F^° m [n] = nf^° m [n\. The functional derivatives (excess-correlation 
potential) of the above relation can be written as [HZl EE] 




(3.40) 



F ex [n(r)} = V H [n(r)} + F xc [n(r)] . 



(3.41) 




(3.42) 



V ex ( 



r, n(r)) 




n— >n(r) • 



(3.43) 
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Information of the homogeneous excess correlation energy /^[n] can be ob- 
tained by subtracting the kinetic energy from the total ground state energy 
of a homogeneous system with N boson. A plot of the DMC data of Magro 
and Ceperley and the STLS data of Strepparola et. al. |115j are depicted 
along with the following fit 

in Fig ()3.6jl . To be consistent with the scaling of Magro and Ceperley [H] 
I will scale all lengths by harmonic oscillator length a = y^H/moo± and the 
energy by e = hu)±/2. Eq. ()3.40|) is solved numerically using the excess 
correlation potentials of Eq. (j3.43J) and Eq. ()3.44|) self-consistently with the 
condition that the areal integral of the condensate density |$(r)| 2 is equal to 
the total number N of bosons. 
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Figure 3.6: Homogeneous excess-correlation energy f ex (p) versus the dimen- 
sionless homogenous density p gas of the VMC data of Magro and Ceperley 
[Hj (pluses) and STLS data of Strepparola et. al. |115] (crosses) compared 
to the fit function Eq. (J3.44)) (dashed-lines). 
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Figure 3.7: Density profiles n{r) for the condensate ( in units a? ho /N ', with 
a ho = ^yh/muj± ) versus the radial distance r (in units Oh ) obtained using 
Eq. iprinjl for iV = 50 and N = 100 bosons 
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4.1 Summary and future directions 



In this thesis I have reported on the development on some fundamental mi- 
croscopic aspects of a two-dimensional Bose-Einstein condensate at finite 
temperature. The approach presented here is of mean-field nature which 
obeys the spontaneous symmetry breaking. This leads to a nonzero mean 
value of the atomic Bose field operator, which corresponds to the order pa- 
rameter for the BEC phase transistion. The order parameter which describes 
the mean-field condensate is highly correlated with quantum fluctuations. 
The condensate which is well described by the Gross-Pitaevskii equation is 
coupled with the Hartree-Fock model of the thermal cloud. We call this the 
two-fluid model description in this thesis. 

Phase fluctuations and scattering properties play a crucial role in the stud- 
ies of 2D Bose gas. One important consequence of lowered dimensionality 
is that the T-matrix for two-body collisions in vacuo at low momenta and 
energy, which should be used to obtain the collisional coupling parameter 
to the lowest order in the particle density, vanishes in the strictly 2D limit 
as the s-wave scattering length becomes larger than the width of the axial 
trapping [101 E2] • It is then necessary to evaluate the scattering processes 
between pairs of Bose particles in by taking into account the presence of a 
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condensate and a thermal cloud through a many-body T-matrix formalism 
IHSl E3j • In fact in this thesis I have calculated the many-body T-matrix 
which includes the effect of thermal excitations. I have found that the many- 
body screening of the interactions due to the occupancy of excited states is 
quite large and rapidly increasing with temperature. However such many- 
body screening has, however, very little effect on equillibrium properties of 
the gas for our choice of system parameters for the temperature of interest in 
this work ( T < 0.5 T c ). This work has been published as "Density profile of 
a strictly two-dimensional Bose gas at finite temperatur", K. K. Rajagopal, 
P. Vignolo and M. P. Tosi, Physica B 344,157 (2004). 

Measurement of Feshbach resonances provides a means of obtaining in- 
formation about interaction strength (repulsive g > or attractive g > 0) 
between atoms. The exact location and width are very important in de- 
ducing the interatomic potentials, but to our knowledge they have not been 
calculated or measured experimentally for the strictly 2D case. Series of cal- 
culation for different species of atomic gasses for the 3D cases can be found 
in the literature |H1 El EH IH31ITU5j. 

In this thesis I have explicitly applied the formalism already developed to 
treat Feshbach resonances in 3D boson gases to the case of a 2D Bose gas 
at zero temperature. One of the main results is that an external magnetic 
field could be used to drive the coupling from repulsive to attractive even in 
a pancake geometry as the scattering length becomes larger than the axial 
width of the condensed cloud. For the case of repulsive interactions I have 
evaluated the particle density profiles at varying coupling strength, demon- 
strating the approach to a Thomas-Fermi profile followed by collapse. In the 
opposite regime of attractive coupling I have given an estimate for the critical 
number of bosons that can form a stable condensate. I found that this critical 
number is inversely proportional to the absolute magnitude of the coupling, 
independently of the strength of the in-plane confinement. This work has 
been published as "Feshbach resonances in a strictly two-dimensional atomic 
Bose gas", K. K. Rajagopal, P. Vignolo and M. P. Tosi, Physica B 353, 59 
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(2004). 

Also in this thesis I have calculated the particle density, ground state 
energy and critical frequency that required to nucleate a single stable vor- 
tex when a 2D condensate put into rotation at various temperatures. My 
calculations illustrate the interplay between the thermal cloud and the struc- 
ture of the condensate in the regions where the radial kinetic energy of the 
condensate is playing an important role (namely, in the outer parts of the 
condensate and in the vortex core) and validated the use of a simple expres- 
sion of the vortex energy for temperatures up to about 0.5 T c . This work has 
been published as " Temperature dependence of the energy of a vortex in 
a two-dimensional Bose gas ", K. K. Rajagopal, P. Vignolo and M. P. Tosi, 
Physics Letter A 328, 500 (2004). 

The work presented in thesis has given a broad microscopic description 
of homogeneous or trapped 2D Bose gas. Thus it opens up an interesting 
range of further research possibilities. At the same time, the genarality of 
the approach considered here leaves numerous questions unanswered. One 
of them which interest me the most is the dynamics of many vortices in a 
2D harmonically trapped Bose gas. In Sec. 3.2 I have developed a density 
functional theory (DFT) technique to calibrate the ground state properties 
of the vortex fluid confined in a two-dimensional harmonic trap within the 
Bogoliubov-de Gennes reference system. 

When a non-interacting auxiliary system is inferred in the DFT model, we 
obtained a Kohn-Sham (KS) type equation. I have calculated numerically 
the density profile and the ground state energy of the vortex fluid using 
KS equation for N < 137 bosons. The excess-correlation potential in this 
equation is obtained from a fit to the liquid phase VMC data of Magro and 
Ceperley jHj and Strepparola et.al |115j for the homogeneous ground state 
energies. 

The crystallization of Yukawa bosons (vortices) can be studied by making 
an analogy to Wigner crystallization |120j in a system of electrons. Wigner 
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first predicted that a system of electrons in a uniform potential would crys- 
tallize at low density within the Hartree-Fock (HF) theory. He proposed 
that at sufficiently low density crystallization greatly reduces the interaction 
energy with only a small increase in the kinetic energy. The energy lowering 
in the crystallization derives from the reduction of the interaction energy by 
spatially separating the electrons. Hartree-Fock studies have shown a tran- 
sition from Fermi liquid to Wigner molecule using small number of electrons 
in a harmonic trap |121l I122j . This transition has been a hot topic in the 
last two years and often associated with the breaking of rotational symmetry 
found in the HF calculation of quantum dots. Very recently Romanovsky 
et. al. |123j have studied crystallization of bosons interacting with Coulomb 
potentials and contact potentials in harmonic trap for small number of par- 
ticles (N = 6). They have demonstrated a two-step method of breaking 
the symmetry of strongly repelling bosons at the HF level followed by a post 
HF symmetry restoration incorporating a correlation beyond mean-field. The 
method describes transition from a BEC state to a crystalline phase in which 
the trap plays acrucial role in localizing the bosons. The Hartree-Fock solu- 
tions breaks the translational invariance of the many-body Hamiltonian and 
bosons are considered to be pinned as in the Wigner crystal. Their method 
can be used to now is to calculate the ground state properties ( density profile 
and energy ) for (N < 137) Yukawa boson in the cystalline phase. 
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